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ABSTRACT 

Simple models of galaxy formation in a cold dark matter universe predict that massive galax- 
ies are surrounded by a hot, quasi-hydrostatic circumgalactic corona of slowly cooling gas, 
predominantly accreted from the intergalactic medium (IGM). This prediction is borne out 
by the recent cosmological hydrodynamical simulations of Crain et al., which reproduce ob- 
served scaling relations between the X-ray and optical properties of nearby disc galaxies. Such 
coronae are metal poor, but observations of the X-ray emitting circumgalactic medium (CGM) 
of local galaxies typically indicate enrichment to near-solar iron abundance, potentially sig- 
nalling a shortcoming in current models of galaxy formation. We show here that, while the hot 
CGM of galaxies formed in the simulations is typically metal poor in a mass-weighted sense, 
its X-ray luminosity-weighted metallicity is often close to solar. This bias arises because the 
soft X-ray emissivity of a typical ~ 0.1 keV corona is dominated by collisionally-excited 
metal ions that are synthesised in stars and recycled into the hot CGM. We find that these 
metals are ejected primarily by stars that form in-situ to the main progenitor of the galaxy, 
rather than in satellites or external galaxies. The enrichment of the hot CGM therefore pro- 
ceeds in an 'inside-out' fashion throughout the assembly of the galaxy: metals are transported 
from the central galaxy by supernova-driven winds and convection over several gigayears, es- 
tablishing a strong negative radial metallicity gradient. Whilst metal ions synthesised by stars 
are necessary to produce the X-ray emissivity that enables the hot CGM of isolated galaxies 
to be detected with current instrumentation, the electrons that collisionally excite them are 
equally important. Since our simulations indicate that the electron density of hot coronae is 
dominated by the metal-poor gas accreted from the IGM, we infer that the hot CGM observed 
via X-ray emission is the outcome of both hierarchical accretion and stellar recycling. 
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1 INTRODUCTION 

Circumgalactic coronae in quasi-hydrostatic equilibrium around 
massive galaxies, formed by the accretion and shock-heating of in- 
tergalactic gas, are a simple but fundamental prediction of galaxy 
formation models in a cold dark matter universe (White & Frenk 
T99T] IKauffmann eTal] [79931 ICole et al.||2000l |Hatton et al 
2003l|Boweret a l. 2006, Croton et al. 2006, Monaco et al.|2007| 



Somerville et al. 2008] |Benson|20 lO). Their properties were first 
explored analytically byjWhite & Fr enk| ( [T991| hereafter WF91), 
who posited that coronae associated with L* galaxies at z = have 
temperatures of T ~ 10 6 K, and emit ~ 10 41 - 10 43 erg s~ x of 
cooling radiation, primarily in the form of soft (~ 0.1 keV) X-ray 
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photons. If present in the local Universe, circumgalactic reservoirs 
of this luminosity should have been detected, if not by previous 
generations of X-ray satellites, certainly by Chandra and XMM- 
Newton. 

The well-documented failure of the Einstein and ROSAT X- 
ray observatories to detect hot atmospheres associated with mas- 
sive disc galaxies (e.g. Bregman & Glassgold 1982; Vogler et al 
[T995] |Bregman & Houck||1997||Fabbiano & Juda||1997||Benson' 
et al. 20001 has often been regarded as a fundamental challenge to 
galaxy formation theory. More recently, the greater sensitivity of 
the Chandra and XMM-Newton observatories has led to detections 
of a hot circumgalactic medium (CGM) associated with both iso- 
lated disc galaxies (Strickland et al. 2004; Wang 2005, Tiillmann 
et al.|2006t|Li et al'|20071 |Jeltema et al.|2008||Owen & Warwick 
2009| [Rasmussen et al. 2009, Sun et al. 2009, L i et al.|20fT) and 
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isolated elliptical galaxies, with no obvious signature of a contribu- 
tion from active galactic nuclei (AGN, David et al. 2006 , Mulchaey 
& Jeltema 2010). However, these data have not established the va- 
lidity of the canonical galaxy formation framework because the in- 
ferred X-ray luminosity is typically one to two orders of magnitude 
fainter than predicted by analytic models. Moreover, the origin of 
the hot gas in these examples remains controversial. 

Most commonly, extra-planar X-ray emission associated with 
disc galaxies is interpreted as a signature of hot ejecta from Type 
II supernovae (SNe) entrained in galaxy- wide outflows. Prima fa- 
cie, this interpretation is not without merit since several studies re- 
port a correlation between the extra-planar X-ray luminosity (Lx) 
and star formation indicators such as Ha line emission or broad- 
band optical/IR luminosity (e.g. Strickland et al. 2004 , Rasmussen 
eT"aLl[2004] |Tullmann et"aT]|2006[ |Sun et al.||2007[ |Li & Wang] 
2012). Particularly spectacular cases such as the archetypal local 



starburst system, M82, bolster this interpretation because of the bi- 
conical morphology of their X-ray surface brightness contours and 
the presence of coincident optical line emission (Strickland et al. 
[2W4|[\^iiieux et al.|2005) . 

In contrast to disc galaxies, the association of hot circumgalac- 
tic gas with elliptical galaxies has been established for many years 
(e.g. |Formanet al.|1979|[T985|[Fabbiano|1989| ). In common with 
disc galaxies, however, the hot gas is commonly interpreted as hav- 
ing an internal origirrjas the ejecta associated with evolved stel- 
lar populations such as Type la SNe and asymptotic giant branch 
(AGB) stars (see e .g.|Mathews|1990||Ciotti et alJig^TJIMathewsl 



|& Brighenti|2003| |Parriott & Bregman|2008^ . This interpretation 
seems reasonable because evolved stellar populations unquestion- 
ably return a significant mass of gas to the interstellar medium 
(ISM) when integrated over the lifetime of the galaxy. Moreover, 
a broad correlation between Lx and the stellar mass of the galaxy 
(AL, or proxies such as Lk) has now been firmly established in 
deep Chandra observations after careful removal of contaminating 
X-ray sources (Boroson et al. 201 1 ). 

The interpretation of diffuse X-ray emission as a signature of 
internal, stellar processes, rather than an accreted coronal reservoir, 
therefore presents a challenge to current models of galaxy forma- 
tion. Motivated by this challenge, we conducted a study (Crain et al. 
2010 hereafter Paper I) using the cosmological hydrodynamic sim- 
ulations of the Galaxies-Intergalactic Medium Interac- 
tion CALCULATION (GIMIC, [Crain et al. 2009. hereafter C09) to 
investigate the origin and nature of the hot CGM associated with 
an unprecedentedly large sample of ~ 500 L* simulated galax- 
ies. Our analysis demonstrated that M82-like systems, whose X- 
ray luminosity is dominated by hot gas entrained in outflows, are 
indeed found in the GIMIC simulations. However, just as in the lo- 
cal Universe, such systems are rare. The majority of ~L* galax- 
ies instead develop extended, quasi-hydrostatic coronae, primarily 
through the shock heating and adiabatic compression of gas ac- 
creted from the intergalactic medium (IGM), supplemented by rel- 
atively small amounts of gas recycled through the galaxy by stellar 
evolution processes. 

This picture appears similar to that posited by the analytic 
model of WF91 but, crucially, the X-ray luminosity of the hot 
CGM of GIMIC galaxies is one to two orders of magnitude fainter 



(at fixed halo mass) than predicted by WF91. The more detailed 
and accurate treatment of gas accretion and feedback in the GIMIC 
simulations highlights that the discrepancy stems primarily from 
inaccurate assumptions regarding the physical structure of the 
corona adopted by WF91. The X-ray luminosity of GIMIC galax- 
ies is therefore broadly consistent with measurements inferred from 
Chandra and XMM-Newton observations of local galaxies. For ex- 
ample, Anderson & Bregman (201 1 ) and Da i et al.|{2012) recently 
reported detections of luminous coronae associated with the giant 
disc galaxies NGC1961 and UGC 12591, respectively, whose in- 
ferred coronal luminosities correspond closely with those of simi- 
larly massive disc galaxies in GIMIC. 

It was noted in Paper I that GIMIC galaxies exhibit a tight cor- 
relation between Lx and halo mass (ALoo), and this naturally gives 
rise to a broad correlation between Lx and both stellar mass (M*) 
and the star formation rate (AL) since, to first order, both of these 
quantities correlate with A/200 • Hence, these observed correlations 
do not necessarily imply an internal origin for the hot gas. This 
raises the appealing possibility of unifying the interpretation of X- 
ray observations of disc and elliptical galaxies in a simple model 
in which, as is inferred in all cosmologically-motivated models of 
galaxy formation, the CGM of galaxies is established by a mixture 
of intergalactic gas accreted throughout their assembly history and 
small amounts of internally recycled gas. 

We intend to explore such a unification in detail elsewhere, 
and focus here on what has long been perceived as a problem for the 
canonical galaxy formation framework: that the hot CGM of ~L* 
elliptical galaxies appears to be enriched to a significant fraction of 
the solar metallicitj ^1 Such high metallicities were first revealed by 
high spectral resolution X-ray spectroscopy of NGC 4636 ( |Xu et al.] 
[2002] ) and NGC 5044 l |Tamura et al.|2003| > using XMM-Newton's 
Reflection Grating Spectrometer (RGS). These measurements have 
been broadly corroborated by imaging spectroscopy with Chandra 
( |Humphrey & Buote|2006] |Athey|2007] hereafter HB06; Ath07), 
which, at the cost of lower spectral resolution, offers a greater ef- 
fective area, wider bandwidth, better isolation of the diffuse gas, 
and the possibility of spatially excising contaminating flux from re- 
solved X-ray point sources. At first sight, the presence of enriched 
hot gas in ellipticals supports a simple 'internal origin' for the gas, 
whereby its metallicity directly reflects that of the evolved stellar 
populations comprising the galaxy (which are approximately near- 
solar, e.g. Trager et al. 2000b a; Terlevich & Forbes 2002 , Gallazzi 
et al. 2006 1, and appears difficult to reconcile with the notion that 
the hot CGM contains a large amount of metal-poor gas accreted 
from the IGM. 

In this paper, we extend the analysis of the hot CGM of z — 
GIMIC galaxies presented in Paper I, to confront the simulations 
with the observational constraints on coronal metallicities presently 
available in the literature. We demonstrate that the simulations are 
compatible with the data, and elucidate the nature of the enrich- 
ment of the hot CGM of > L* galaxies, demonstrating how, where 
and when the metals were injected. Building on the arguments pre- 
sented in Paper I, we consider this compatibility to be further ev- 
idence that X-ray data support the canonical view of galaxy for- 
mation in the CDM paradigm, in a fashion that is orthogonal to, 
for example, the confrontation of galaxy formation models with 



1 Since we are mainly concerned here with ~L* galaxies, we do not con- 
sider galaxies with significant AGN activity, as were considered theoret- 
ically by e.g. |Ciotti & Ostriker| j!997| [2007}, nor those confined by the 
potential of galaxy groups. 



2 The first metallicity measurement of the hot CGM of a late-type galaxy 
(NGC891) was reported by |Hodges-Kluck & Bregman| {2012) during the 
preparation of this manuscript. They report Z ~ 0.1 Z©, consistent with 
the gas being primarily sourced by accretion from the IGM. 
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optical/infra-red data (e.g. the stellar mass function or the Tully - 
Fisher relation). This paper is laid out as follows. In §[2]we provide 
a brief overview of the GIMIC simulations and our analysis of them. 
We confront these analyses with X-ray data in §15] and demonstrate 
how the enrichment proceeds in § H] Finally, we discuss and sum- 
marise our findings in § B] 



2 METHODS 

2.1 Simulations 

The GIMIC simulations are described in detail in C09, where thor- 
ough discussions of the generation of the initial conditions, the sim- 
ulation algorithms, and initial results may be found. Techniques for 
computing gas-phase X-ray luminosities and stellar optical lumi- 
nosities are described in detail in Paper I. Here, we present only a 
brief overview and limit the description to aspects that are specifi- 
cally relevant to this study, referring the reader to C09 and Paper I 
for additional details. 

Using 'zoomed' initial conditions, GIMIC follows with full gas 
dynamics and, at relatively high resolution for cosmological vol- 
umes, the evolution of five roughly spherical regions of comov- 
ing ra dius r ~ 20 h~ M pc, drawn from the Millennium Simula- 
tion ( |Springel et al.|2005| . The remainder of the 500 3 (/i _1 Mpc) 3 
Millennium Simulation volume is modelled with collisionless par- 
ticles at much lower resolution, thus correctly following the cosmic 
large-scale structure. In order to trace a wide range of cosmic en- 
vironments, the regions were chosen such that their overdensities 
deviate by (—2, — 1, 0, + 1,+2)<t from the cosmic mean, where a 
is the rms mass fluctuation on the 20 hT 1 Mpc radial scale of the 
regions, at z = 1.5. 

The GIMIC initial conditions were generated at intermedi- 
ate resolutiorr] (m gas = 1.16 x 10 7 ft~ 1 M Q ) and high res- 
olution (m gas = 1.45 x 10 6 hT 1 Mo), and ran with gravita- 
tional force softening for the high-resolution particles of e£wj = 
(1.0, 0.5) h~ kpc for the intermediate- and high-resolution runs 
respectively. The relatively large cosmological volume traced with 
gas dynamics produces a sample of over 600 relatively isolated 
~L* galaxies, making GIMIC complementary to ultra-high reso- 
lution simulations (m gas ~ 10 4 hT 1 Mq) of individual galaxies. 

The cosmological parameters adopted for GIMIC are the same 
as those assumed by the Millennium Simulation and correspond to 
a ACDM cosmogony with fi m = 0.25, tt A = 0.75, O b = 0.045, 
a s = 0.9, Ho = lOO/ikms^Mpc- 1 , h = 0.73, n s = 1 (where 
n s is the spectral index of the primordial power spectrum). The 
value of (78 is approximately 2<r higher than inferred from the most 
recent cosmic microwave background (CMB) and large-scale struc- 
ture data (Komat su et al.|20TT) , which will affect the abundance of 
Li, galaxies somewhat, but should not significantly affect their in- 
dividual properties. 

The initial conditions were evolved from z — 127 to z = 
with the TreePM-SPH code GADGET3, a substantial upgrade of 
GADGET2 (Springel 2005) that features improved load-balancing 
2008 ) and incorporates a wide array of new 



(as in Springel et al 
baryon physics modules. Radiative cooling rates are computed on 
an element-by-element basis by interpolating pre-computed tables, 
generated with CLOUDY (version 07.02, |Ferland et al.|1998| >, that 



3 We reserve the term 'low resolution' for the original realisation of the 
Millennium Simulation. 



specify cooling rates as a function of density, temperature and red- 
shift, and that account for the presence of the CMB and a spatially- 
uniform Haardt & Madau (2001) ionising UV/X-ray background 
(for further details, see Wiersma et al. 2009a). The ionising back- 
ground is imposed at z = 9 and effectively quenches star forma- 
tion in small haloes (Oka moto et al.| [2008 C09). Star formation 
is implemented following the prescription of Schaye & Dalla Vec- 
[chia|J2008> 



(Kennicutt 



which reproduces the observed Kennicutt-Schmidt law 



1998) by construction. 



Feedback from SNe is implemented using the kinetic wind 
model of Dalla Vecchia & Schaye ( 2008 1, with an initial wind ve- 
locity v w = 600 km s~ and mass-loading parameter (the ratio 
of mass of gas that receives an impulsive kick to the mass of stars 
formed) of /3 = 4. This scheme is energetically feasible, requiring 
approximately 80 percent of the total energy available from Type 
II SNe for our adopted choice of initial mass function (IMF), that 
due to Chabrier (2003). This choice of parameters produces a good 
match to the peak of the cosmic star formation history (see C09, 
Schaye et al. 2010). The timed release of individual elements by 
massive stars (Type II SNe and asymptotic giant branch stars) and 
intermediate-mass stars (Type la SNe) is incorporated following the 
prescription of Wiersma et al. (2009b). A set of 1 1 elements (H, He, 
C, Ca, Fe, Mg, N, Ne, O, S, Si), representing the most important ra- 
diative coolants, are followed, enabling us to determine the element 
abundances of individual SPH and star particles. 

The GIMIC simulations are particularly well suited to the study 
of ~L* galaxies. As shown in Paper I, the implementation of ef- 
ficient (but energetically feasible) feedback from SNe largely pre- 
vents overcooling on the mass scale of L* galaxies, and is key to 
the reproduction of the observed X-ray scaling relation presented in 
that study. Indeed, GIMIC accurately reproduces the rotation speeds 
and star formation efficiencies of z = disc galaxies for 10 9 < 
M* < 10 10 - 5 M , although galaxies with M* > 10 11 M Q do still 
suffer from some overcooling (McCarthy et al. 2012b). Moreover, 
|Font et aL] J20TT} demonstrated that L* galaxies in GIMIC exhibit 
satellite luminosity functions and stellar spheroid surface bright- 
ness distributions that are comparable to those of the Milky Way 
and M31, whilst McCarthy et al. (2012a) further demonstrated that 
this correspondence extends also to their global structure and kine- 
matics. 



2.2 Galaxy & halo identification 

We identify bound haloes using the SUBFIND algorithm presented 
by Dolag et al. (2009), which extends the standard implementa- 
tion (Springel e t al.||206T| > by including baryonic particles when 
identifying self-bound substructures. This procedure first finds dark 
matter haloes using a friends-of- friends (FoF; |Davis et al.| (jl985)) 
algorithm, with the standard linking length in units of the inter- 
particle separation (6 = 0.2) It also associates all baryonic parti- 
cles with their nearest neighbour dark matter particle. The aggre- 
gated properties of baryonic particles associated with grouped dark 
matter particles define the baryonic properties of each FoF group. 
Halo substructures are then identified using a topological unbind- 
ing algorithm. This provides an unambiguous definition of a galaxy 
within the simulations, namely the set of star particles bound to in- 
dividual subhaloes. The gas bound to each subhalo then forms the 
ISM and CGM; these phases are delineated by the density thresh- 
old for the onset of star formation, nu = 0.1 cm" . Since individ- 
ual haloes may have more than one 'subhalo', they may host more 
than one galaxy. In analogy to semi-analytic models, the stars as- 
sociated with the most massive subhalo of a FoF group are defined 
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as central galaxies, whilst stars associated with substructures are 
satellites. When referring to the mass of haloes, we adopt spherical 
overdensity (SO, Lacey & Cole|1994) l masses, centred on the most 
bound particle within FoF groups, with radius such that the mean 
enclosed density is 200 times the critical density, i.e. M200. 



2.3 Simulated galaxy sample 

As in Paper I, we draw galaxies from all five intermediate- 
resolution GIMIC regions, and use the high-resolution — 2<r region 
as a benchmark for brief convergence tests, which are presented 
in Appendix [A| We impose a number of selection criteria in order 
to construct a sample similar to observational samples: we select 
central (i.e. the most massive galaxy within an FoF halo) galaxies 
whose stellar mass lies in the range 10 10 < KL < 10 115 M Q , 
corresponding to a fairly broad mass range bracketing L* . We ad- 
ditionally require that the subhalo with which each galaxy is as- 
sociated should account for at least 90 per cent of the mass of the 
parent FoF group, i.e. M su b > 0.9Mf o p. This criterion is intended 
to select only isolated systems by excluding those that are interact- 
ing or are members of galaxy groups or clusters. These selection 
criteria result in a sample of 617 galaxies at z = 0. 

As shown in |McCarthy et al.| ( |2012b) , the intermediate- 
resolution GIMIC simulations accurately reproduce the stellar mass 
- rotation speed (or 'Tully-Fisher') relation at z — for galax- 
ies with 10 9 < M* < 10 10 ' 5 Mq. These galaxies are consistent 
with the M*-M2oo relation inferred from stacked weak lensing + 
satellite kinematics analyses (see Dutton et al.|2010| and references 
therein). [McCarthy et al.| ( |2012b| also found that the most massive 
galaxies in our sample (M* > 10 11 Mq, corresponding to A'-band 
luminosities Lk > 10 11 Lk,©) exhibit some overcooling, since our 
feedback scheme becomes inefficient within massive haloes (see 
C09). We retain these galaxies in our sample in order to connect 
with observations of similarly massive ellipticals, but caution that 
these galaxies are likely hosted by dark matter haloes that are in- 
sufficiently massive, and exhibit star formation rates at low-redshift 
that are too high. 



2.4 Computing optical & X-ray luminosities 

Optical and X-ray luminosities are computed in postprocessing. 
The former are calculated individually for star particles, consid- 
ering them as simple stellar populations (SSP). Star particles are 
formed with the IMF due to Chabrier (2003), and their spec- 
tral energy distribution (SED) is derived by interpolation over the 
GALAXEV models of |Bruzual & Char iot (2003). The optical lu- 
minosity of each particle is obtained by integrating the product of 
its SED with the appropriate filter transmission function; the over- 
all broadband luminosity of a galaxy is then defined as the sum of 
the luminosities of all star particles comprising a galaxy. Since we 
are primarily concerned here with if-band luminosities, we do not 
expect dust extinction to be an important consideration. 

Gas phase X-ray luminosities are also computed on a per- 
particle basis, following: 



Lx — n c riiAV, 



(1) 



where n c and ni are the number densities of electrons and ions re- 
spectively, A is the cooling function in units of erg s" 1 cm" 2 (inte- 
grated over an appropriate passband, for example 0.5-2.0 keV), and 



Table 1. Abundances of the 1 1 elements tracked self-consistently within 
the GIMIC simulations, and used to compute the radiative cooling rate and 
X-ray luminosity of gas, presented for the Anders & Grevesse 1 19891 abun- 
dance standard (assumed by the APEC model), and that adopted by the 
GIMIC simulations, which is derived from CLOUDY (version 07.02). We 
follow the convention log 10 Nn = 12. 



Element 


AG89 


Cloudy 


Hydrogen 


12.00 


12.00 


Helium 


10.99 


11.00 


Carbon 


8.56 


8.39 


Calcium 


6.36 


6.36 


Iron 


7.67 


7.45 


Magnesium 


7.58 


7.54 


Nitrogen 


8.05 


7.93 


Neon 


8.09 


8.00 


Oxygen 


8.93 


8.69 


Sulphur 


7.21 


7.26 


Silicon 


7.55 


7.54 



table generated using the Astrophysical Plasma Emission Code 
(APEC, vl.3.1, see |Smith et al.p OOT) under the assumption that 
the gas is an optically thin plasma in collisional ionisation equilib- 
rium. APEC cooling rates are computed on an element-by-element 
basis for each particle, 



A J (T J ) = -^^^A fc (T J ), 



(2) 



where Tj is the gas temperature of the j th particle, Zj^ is the metal 
mass fraction of the element k for this particle, and \k(Tj) is the 
cooling function for the element k. The total cooling rate is arrived 
at by summing over the TV = 1 1 most important elements for cool- 
ing (H, He, C, Ca, Fe, Mg, N, Ne, O, S, Si), which are individually 
and self-consistently tracked during the simulation. APEC assumes 
the solar abundance ratios of Anders & Grevesse ( 1989. hereafter 
AG89), but it is straightforward to rescale the cooling rates for 
arbitrary abundances, which we do here on a particle-by-particle 
basis, normalising by the metal mass fractions assumed by AG89 
(Zk^AGsg), to achieve consistency with our chemodynamics im- 
plementation. For completeness, we present the solar abundances 
adopted by AG89 and GIMIC in Table [T] The latter corresponds 
to the abundances assumed by CLOUDY (version 07.02, last de- 
scribed by Ferland et al. 1998). Our scheme excludes gas within 
the ISM (particles with non-zero star formation rate, see Schaye & 
Dalla Vecchia 2008) since we assign these particles a temperature 
ofT= 10 4 K. 



2.5 Computing iron abundances 

Since iron has the strongest diagnostic lines in the soft X-ray band, 
Zp B is the most commonly used proxy for gas-phase metallicities 
in low temperature plasmas (Loewenstein & Mushotzky 19981. We 
therefore adopt the mass-weighted iron mass fraction as a metallic- 
ity proxy, computed thus: 



ymw 



£f^ 



FejTOj 



£ 



(3) 



V = m gas /p, where m gas is the gas particle mass and p the gas where Zf c j is the SPH kernel-weighted iron abundance of the j 



particle density. We compute A by interpolating a pre-computed 



gas particle (see discussion in §12. 6b, and rrij is its mass. The sum 
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runs over the N gas particles associated with each galaxy that have 
a non-zero soft X-ray luminosity (and thus comprise the hot CGM, 
see § [24). 

In Paper I, we cautioned that the sensitivity of the X-ray emis- 
sivity of gas to density, temperature and metallicity can lead to 
an unrepresentative view of the overall state of hot gas. In more 
specific terms, the unavoidable luminosity-weighting of X-ray ob- 
servations can yield measurements that differ markedly from an 
ideal, mass-weighted measurement. The most direct comparison 
between simulated and observed coronae is obtained by deriv- 
ing Lx-weighted quantities in the former. Therefore, we also de- 
fine a luminosity-weighted metallicity indicator, that is appropriate 
when comparing with observationally-inferred metallicities, by ex- 
changing the particle masses (rrij) in Eqn.[3]for their SPH kernel 
weighted soft X-ray luminosities (Lxj): 



„i w _ Ef Zfc, 3 L x ,j 

Is] L *,j 



^Fc 



(4) 



In general, we use the mass weighted measure to assess the physi- 
cal state of the hot CGM, and reserve the luminosity weighted mea- 
sure for direct comparsions with observationally-inferred measure- 
ments. 



2.6 Metal mixing 

The cooling rate, and therefore star formation rate, of cosmic gas 
is strongly affected by its metallicity, so the mixing of stellar nu- 
cleosynthetic products into the ISM and CGM is clearly of signif- 
icance for models of galaxy evolution. The process of metal mix- 
ing, however, remains poorly described by theory, and poorly con- 
strained by observation. 

Nonetheless, a common theme amongst both types of study 
is that mixing is not particularly efficient, across a broad range of 
scales. On relatively small scales (a few parsecs), the Eulerian sim- 
ulations of Parriott & Bregman (2008 1, and particularly Bregman & 
Parriott ( 2009 ), indictate that much of the ejecta from evolved stars 
and planetary nebulae remain unmixed with the hot ISM/CGM. 
This picture appears to be corroborated by observations of the 
structure of stellar ejecta at both infra-red (Wareing et al.||2006) 
and ultra-violet (Martin et al. 20071 wavelengths. On significantly 
larger scales, Schaye et al. ( 2007 1 argue that the large population of 
compact, intergalactic CIV absorbers seen in their sample of QSO- 
sightlines is evidence for inefficient mixing of metals entrained in 
winds and transported from the ISM into the CGM/IGM. 

The physical interpretation of the distribution and transport of 
metals seen in our simulations therefore requires an understanding 
of the behaviour and limitations of the numerical implementation 
of these processes. The formulation of SPH naturally defines three 
spatial scales for which a discussion of metal mixing is appropriate: 
scales smaller than, comparable to, or larger than the SPH kernel 
(i.e. the smoothing length). 

On the smallest scales, the nucleosynthetic products of stellar 
evolution are implicitly assumed to be completely mixed, since the 
metal mass released by the SSP represented by each star particle is 
assigned to neighboring SPH particles. Clearly, each particle is the 
smallest unit over which mass can be distributed, so it is therefore 
possible that on subgrid scales mixing is overestimated. Such scales 
are, by definition, below the resolution limit of the simulation, and 
as such physical insight should not be drawn from the simulations 
on these scales. 

On scales comparable to the kernel, the finite numerical sam- 



pling of the fluid, and the absence of implicit diffusion of scalar 
quantities associated with SPH particles, suppresses the mixing of 
entropy and metal^] (Wad sley et al.| |2008). This shortcoming is 
clearly of relevance here, and it has been shown that the metal- 
licity of circumgalactic gas in simulations is increased when sub- 
grid models for efficient turbulent mixing are imposed ( She n et ak] 
2010). However, we opt against imposing such a scheme here, since 
it represents represent an attempt to resolve a numerical problem 
with a model of a poorly understood physical process. 

We intend to explore small scale diffusion in future work, 
using SPH schemes augmented with higher-order dissipation 
switches (e.g. Read & Hayfield 2012) that diffuse scalar quantities 
implicitly. Here, we follow Wiersma et al. (2009b) and minimise 
the sampling problem by computing cooling rates, X-ray luminosi- 
ties and overall coronal iron mass fractions (Eqns. [3] & HI using 
kernel-weighted element abundances averaged over a given SPH 
particle's neighbours. The resulting smoothed metallicity comple- 
ments the 'particle metallicity' that is the sum of the element mass 
fractions carried by each particle. This process has the potential to 
bias the recovered metallicity of hot circumgalactic particles, be- 
cause smoothed metallicities are computed using all SPH neigh- 
bour particles, irrespective of their thermal state. The smoothed el- 
ement abundances of hot particles close to the disc-corona interface 
are therefore computed from both cold (i.e. ISM) and hot (coro- 
nal) particles. We gauge the degree to which mixing on the scale 
of the kernel influences our conclusions, by comparing results de- 
rived from particle and smoothed metallicities. We conclude from 
this analysis, presented in AppendixfB] that mixing on the scale of 
the kernel has no bearing on our findings. 

On larger scales still, simulations indicate that the distribution 
of metals is governed chiefly by winds driven by energetic feed- 
back (jWiersma et al. 2011), a process self-consistently modelled 
by our simulations. Whilst the standard implementation of SPH is 
well-understood to suppress the development of Rayleigh-Taylor 
and Kelvin-Helmholtz instabilities that would assist the mixing of 
metals entrained in winds into their ambient medium, the inefficient 
mixing inferred by Schaye et al. (2007) on these relatively large 
scales suggests that it is reasonable to conclude that predictions of 
the large-scale distribution of metals derived from our simulations 
are robust to the implementation of metal mixing adopted by our 
simulations. It is the metal distribution on these large scales that is 
the central the focus of this study. 



2.7 Stacked maps 

To illustrate the correspondence of the simulated galaxies with 
observable properties, we present in Fig. [T] maps derived from 
stacking the 148 galaxies from our sample whose halo masses 



lie in the range 10 



< M 2O o[M ] < 10 1 



The maps are 



500 kpc across, with the exception of the bottom left panel, which 
is zoomed-in to 50 kpc, and in each case the galaxies have been 
rotated into the edge-on orientation by aligning the vertical axis 
with the angular momentum vector of all stars within 35 kpc of 
the galactic centre. The left hand column shows the soft X-ray sur- 
face brightness (E0.5-2.0kev) of the hot CGM, the central column 
shows the surface density of the hot CGM (Ehot-CGM; top) and 
of stars (E,; bottom), whilst the right hand column shows the iron 



4 It is not clear that Eulerian algorithms present a more appealing alterna- 
tive, since they can in fact mix gas too efficiently (e. g. [Mitchell et al. 2009). 
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Figure 1, Stacked maps of the baryons associated with the 148 GIMIC galaxies (from our sample of 617) with halo mass in the range 10 12 - 25 < Af2oo/MQ < 
10 12 ' 50 . The galaxies have been rotated so that the angular momentum of stars is parallel to the vertical axis prior to stacking. The maps are 500 kpc across, 
with the exception of the bottom left panel, which is zoomed to 50 kpc. Left column: the soft X-ray (0.5-2.0 keV) surface brightness. Centre column: the 
surface density of hot gas (top) and stars (bottom). Right column: the iron abundance of hot gas (top) and stars (bottom). The X-ray emission is markedly 
more centrally peaked than the mass distribution of the hot CGM, owing to its strong dependence on density, metallicity and temperature. A halo-wide iron 
abundance gradient is evident in the hot CGM, that does not trace the underlying distribution of stars, indicating the prevalence of metal transport. 



abundance of the hot CGM (^Fo,hot~CGM; top) and of stars (Zp e ,*; 
bottom). 

The alignment procedure is clearly effective: a disc can be 
seen in the stellar surface density map, as can the spheroidal stellar 
halo component that was explored by Font et al.lpOTTb and Mc- 
|Carthy et al.| ( |2012a) . Similarly, the stellar iron abundance shows a 
clear disc plane, and a large-scale gradient similar to the metallicity 
gradients observed in nearby stellar haloes (Foster et al. 2009 , Spo- 
|laor et al.|2010| >. The presence of strong disc features in these maps 
is unsurprising, as roughly two-thirds of the galaxies comprising 
our sample are disc-dominated (see Paper I). The hot gas maps ex- 
hibit circular surface density and soft X-ray surface brightness dis- 
tributions, indicative of a quasi-spherical hot CGM associated with 
the galaxies. The coronal iron distribution is centrally peaked, but 
does not trace the underlying stellar iron abundance, indicating that 
metal transport is important for enriching the hot CGM. We return 
to the issue of metal transport in §E] 

A qualitative difference in the radial distribution of the coro- 
nal gas and its X-ray luminosity is evident in the maps: the former 
is much more extended than the latter. This follows from the com- 
plex mapping of the state of the ionised plasmas (p, T, Z) into an 
X-ray emissivity: bremsstrahlung radiation is weakly dependent on 
temperature (oc T 1 ' 2 ) and strongly dependent on density (oc p 2 ), 



whilst the metal line emission that dominates the emissivity of low 
temperature plasmas is also dependent on the abundance of metals. 
A quantitative assessment of the difference is presented in Fig. [2] 
where we bin the stacked maps radially, and plot the radial surface 
density of the hot CGM (black) and its X-ray surface brightness 
(red). Solid lines show the median value of each quantity, arbitrar- 
ily normalised to peak at a value of unity, whilst shaded regions rep- 
resent the scatter between the 10 th and 90 th percentiles. The inset 
plot presents the same data on logarithmic axes, and the overplot- 
ted dashed green lines represent fits, using the commonly adopted 
/^-function: 

(0.5-3/3) 



E(r) oc 



(5) 



to the stacked surface density and surface brightness profiles, where 
r c is the core radius and /? characterises the profile shape. It is 
immediately evident that X-ray emission does not trace the hot 
gas mass in a simple fashion: the surface brightness profile is 
much more centrally peaked than the surface density profile. This 
is reflected quantitatively in the least-squares reduced parameters 
of the fitted /^-functions: the surface density is best described by 
(r c , f3 = 120 kpc, 0.59), whilst the surface brightness is best de- 
scribed by (r c ,/3 — 12 kpc, 0.56). The recovered f3 values are 
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Figure 2. Circularly-averaged median projected radial profiles of the sur- 
face density of the hot CGM (black) and its soft X-ray surface brightness 
(red), derived by radially binning the maps presented in Fig. [T] and nor- 
malising to the maximum value of the median in each case. Shaded regions 
represent the scatter between 10 th and 90 th percentiles, and the inset plot 
presents the same data on logarithmic axes. Green dashed lines represent 
best fit /3 functions to both profiles (see text for parameters). It is clear that 
the surface brightness is much more centrally peaked than the surface den- 
sity, and therefore that X-ray emission does not trace the hot gas mass in a 
simple fashion. 



consistent with those derived from a Monte Carlo Markov Chain 
analysis of stacked ROSAT pointings of 2165 galaxies from the 
2MASS Very Isolated Galaxy Catalogue, recently presented by An- 
|derson et al. 1(20 12) , and similar to the values derived from deep 
observations of NGC 1961 and UGC 12591 by |Anderson & Breg-| 
|man| pOTT) and |Dai et al.| (2012), respectively. Those authors in- 
fer 1 kpc cores for NGC 1961 and UGC 12591, and constrain the 
stacked core radii to 0.5 — 5.0 kpc. This is somewhat smaller than 
the core radius we recover from our stacked surface brightness pro- 
file, highlighting a potential shortcoming of the simulations and an 
interesting avenue for future study. However, having also fitted /?- 
functions to the individual profiles of the well-sampled galaxies in 
our sample, find that a wide range of core radii are recovered; a 
small but significant fraction of the profiles (~ 15 percent) are in 
fact best described by cusped profiles (r c — > 0). The simulations 
therefore yield a non-negligible fraction of galaxies that are similar 
to those observed. Moreover, we caution that the small sample of 
observationally inferred fits may be unrepresentative of the popula- 
tion we explore with the simulations, since galaxies with small core 
radius exhibit a greater peak surface brightness at fixed luminosity, 
rendering them easier to detect (c.f Rasmussen et al. 2006). 

Besides highlighting that it is challenging to observe hot coro- 
nae at radii significantly beyond the optical radius of the central 
galaxy, Fig. [2] indicates that, in the absence of spatially-resolved 
spectroscopy (which is extremely challenging to realise for the hot 
CGM of galaxies with current instrumentation), the intrinsic lu- 
minosity weighting of observations biases the coronal metallicity 
inferred from X-ray spectroscopy towards the high surface bright- 
ness, metal rich gas close to galactic centres. The importance of this 
effect will be explored in the following section. 



3 THE CORONAL METALLICITIES OF PRESENT-DAY 
GALAXIES 

In this section, we present the iron mass fractions of the hot CGM 
associated with GIMIC galaxies at z = 0, and confront them with 
observational measurements inferred from X-ray imaging spec- 
troscopy. We start with a brief discussion of observational studies 
in the literature. 



3.1 Confrontation with observational data 

The excellent spatial resolution of Chandra has revolutionised 
studies of metal enrichment in galaxies and galaxy groups. The 
so-called 'Fe-discrepancy' (Arimoto et al. 1997, Renzini 1997), 
namely the inference of highly subsolar iron abundances in early 
type galaxies, has been shown to result from the inability of ROSAT 
and ASCA to resolve coronal temperature gradients (HB06), that 
can lead to underestimated abundances when isothermal spectral 
fits are adopted (e.g. Buote & Canizares 1994, Trinchieri et al. 



|1994| |Buote & Fabian||1998[ l. Early studies with Chandra (and 
XMM-Newton) focussing on individual systems indicated near- 
solar coronal iron mass abundances (e.g. |Buote|2 002, Gastaldello 
|& Molendi| [20021 |Buote et all [2053 ; O'Sulli vanet al.||20031 
Humphrey et al. 2004, Kim & Fabbiano 2004). Motivated to ex- 

I 1 ! II I I II I 

tend this discovery to a representative sample of galaxies, HB06 
uniformly reduced archival Chandra data for 28 nearby, early type 
galaxies with soft X-ray luminosities spanning the range 10 39 < 
Lx [erg s _1 ] < 10 43 . They report that their measurements support 
an absence of correlation between Z-p c and Lx, and confirm that in 
all but one case (NGC 1553), the coronal iron abundance corre- 
sponds closely with that of the stars. A similar study was presented 
by Ath07, who uniformly reduced Chandra imaging spectroscopy 
data (from the ACIS-S3 instrument) for 54 nearby (d < 110 Mpc) 
early type galaxies. For the brightest 22 galaxies in this sample, 
Ath07 inferred iron abundances via spectral modelling adopting 
variable abundance ratios, and found that only a narrow spread of 
ratios were required for good fits. The lower signal-to-noise ratio of 
the remaining 32 fainter galaxies required that their metallicities be 
derived subject to the assumption that their abundance ratios follow 
those of the brighter galaxies. In common with HB06, Ath07 report 
typical iron abundances that are approximately solar. We note that 
22 galaxies are present in both samples. 

In this study we adopt the HB06 sample of 28 galaxies, and 
supplement it with the 32 galaxies from Ath07 that are not present 
in HB06; in general these are fainter systems than were studied by 
HB06. For reference, we show briefly in Appendix [C] how closely 
the Zfc and Lx measurements inferred by these two studies match. 
Whilst we refer the reader to HB06 and Ath07 for detailed descrip- 
tions of the data reduction, it is appropriate to discuss key aspects 
of their analyses. Firstly, Ath07 adopt the superceded solar abun- 
dance standard of Anders & Grevesse ( 1989), whilst HB06 adopt 
the more recent standard of Asplund et al. (2005), which is very 
similar to that adopted by GIMIC. For consistency, in what follows, 
we therefore multiply the iron abundances derived by Ath07 by 
a factor of 1.65 to achieve consistency between the adopted stan- 
dards. Where the signal-to-noise ratio of the X-ray spectroscopy 
allows, both studies fit plasma models with variable abundance ra- 
tios; in the event of insufficient S/N, HB06 fix to solar ratios, whilst 
Ath07 adopt the typical ratios of the high-S/N galaxies. HB06 re- 
port estimated 'global abundances', whereby systems with mea- 
sureable abundance gradients are fit by a power-law and extrap- 
olated to 30 arcminutes. Ath07 report hot gas abundances within 
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Figure 3. The X-ray luminosity-weighted iron abundance of circumgalactic gas, as a function of soft X-ray luminosity (left) and the i\~-band luminosity of 
galaxies (right). The left-hand panel comprises the 28 nearby elliptical galaxies analysed by HB06 (circles) and the 26 (from a total sample of 54) nearby 
ellipticals from the sample of Ath07 (triangles) that do not overlap with the HB06 sample. The right hand panel comprises the subset of the overall galaxy 
sample for which K-bund magnitudes are available in the 2MASS Large Galaxy Atlas. Open circles with upward arrows represent cases where the HB06 
spectral analysis yielded no formal upper limit on the iron abundance. Symbols are overplotted on the median Lx-weighted Zp c — Lx an d ^Fc — ^K 
relations of galaxies in the GIMIC simulations, with la and 2<r scatter shown by the densely- and sparsely-hatched red regions respectively. Although the 
correlation between Zp c and Lx, implied by the addition of the Ath07 sample to the HB06 sample, is not reproduced by the simulation (but see Appendixlcl, 
the normalisation of the observed metallicities is broadly reproduced by the simulations. 



one effective (optical) radius of the host elliptical. Finally, we note 
that the spectral analysis performed by HB06 produced no upper 
bound on the iron abundance of 12 of 28 galaxies in their sample; 
we plot these are open circles with no upper error bar. 

We compare the Zfc-Lx relation of the HB06 and Ath07 
samples with the GIMIC results in the left-hand panel of Fig. [3] 
The observational measurements are overplotted on the median 
and la and 2a scatter (represented by the densely- and sparsely- 
hatched regions, respectively) of the Lx-weighted coronal iron 
abundance of the 617 galaxies in the GIMIC ~L* galaxy sample. 
The GIMIC sample does not extend to such high soft X-ray lumi- 
nosities (> 10 42 erg s _1 ) as the data, since overcooling in haloes 
with these luminosities produces galaxies that are much more mas- 
sive than Li, at z = and are hence excluded from our sample. 

The simulated galaxies exhibit no strong correlation between 
the Lx-weighted iron abundance and Lx, with a median Zp c that 
is near-solar across five decades in Lx, and a 2a-scatter of approx- 
imately one decade in Zp c . This absence of correlation suggests 
that simulations with efficient feedback in massive haloes would 
not generate significantly different iron abundances to those shown 
here. In fact, McCarthy et al. (2010) showed that including efficient 
feedback in galaxy groups, implemented as AGN feedback in the 
OWLS simulations (Schaye et al. 20101, has almost no impact upon 
the iron injected into the hot intragroup/intracluster medium, even 
though the mass of stars in galaxy groups/clusters is supressed by a 
factor of four. This is because the additional iron that is synthesised 
in simulations without AGN feedback is mostly injected into the 
cold ISM (rather than the hot CGM), and is rapidly reincorporated 
into subsequent generations of stars. 

HB06 concluded that their sample shows no convincing evi- 
dence of a correlation between Zp c and Lx- The addition of the 
Ath07 data, which populates the Zf c -Lx plane at lower luminosi- 
ties, indicates an increase in metallicity from low to high lumi- 



nosity systems. The simulated sample does not exhibit such a pro- 
nounced correlation. Whilst this potentially represents a shortcom- 
ing of the simulations, we note that the correlation remains ill- 
understood observationally and, as discussed by Buote & Fabian 
( 1998 ) and HB06, the recovery of Zp c is potentially hampered by 
complex, multi-component thermal structure in the hot gas. This 
tends to bias the measurement towards low abundances. A second 
potential shortcoming of the model is the absence of galaxies with 
Zp,, > 2 Zq, as are present in the HB06 sample. We again caution, 
however, that the quoted uncertainties potentially underestimate the 
true uncertainy of each measurement. As an example, we note that 
NGC 4365 features in both samples: HB06 infer a metallicity of 
Zp c — 5 Zq, whilst Ath07 assigns it a sub-solar metallicity. This 
uncertainty is discussed further in Appendix |C| 

We stress that it is not our aim to advance the GIMIC sim- 
ulations as a complete and detailed description of the CGM as 
observed in soft X-rays. Whilst a detailed dissection of the estab- 
lishment of the Zp c - Lx relation is clearly an interesting avenue 
for future study, it is our aim here to test whether observational 
measurements rule out the general framework within which the 
hot CGM is established primarily via the accretion and heating of 
metal-poor intergalactic gas, as proposed in Paper I, The left-hand 
panel of Fig. [^demonstrates conclusively that the presence of a hot 
CGM, dominated by low-metallicity gas accreted from the IGM, 
is perfectly consistent with luminosity-weighted, near-solar metal- 
licities. In our cosmological simulations, in which most of the hot 
CGM is of extragalactic origin, the median luminosity-weighted 
iron abundance of the gas is typically ~ 0.4 — 0.7 Zq, with 2a 
scatter between ~ 0.2 — 1.1 Zq, across five decades in X-ray lumi- 
nosity. Thus, our simulations demonstrate that the often-used argu- 
ment that high metallicities preclude a predominantly extragalactic 
origin for the hot CGM of typical galaxies is not robust. 

As a consistency check we plot, in the right-hand panel of 
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Fig. [3] the Zfb-Lk relation for the subset of galaxies in the HB06 
and Ath07 samples for which if -band magnitudes ar e available in 



10.0 



the 2MASS Large Galaxy Atlas ^Jarrett et al.|2003| 



The data ex- 
hibit a larger scatter (approximately two decades in Zp c at Lk— 
10 Lk,©), which is broadly consistent with the simulations. The 
simulations indicate a weak correlation that is broadly consistent 
with the observational measurements over this luminosity range. 
The Lx-Lk relation for GIMIC galaxies was presented in Paper I, 
where it was shown to accurately reproduce the observed relation 
for disc galaxies. We intend to present a comparison with both disc 
and elliptical galaxies in a forthcoming paper (Crain et al. in prep). 
There is a reasonable match between the results of the sim- 
ulated galaxies and the observational measurements: nearly all of 
the measurements lie within the 2<r scatter of the simulation re- 
sults, over five decades in Lx- This is remarkable given that the pa- 
rameters of the subgrid modules in the simulations were not tuned 
to reproduce any particular X-ray scaling relations. We conclude 
that models of galaxy formation based upon the CDM paradigm in 
which the hot CGM of > L* galaxies is dominated by gas accreted 
from the IGM, are consistent with the best available observational 
constraints on the heavy element abundance of hot gas associated 
with such galaxies. 



3.2 Recycled mass and luminosity profiles 

As motivated in the Introduction, the correspondence of near-solar 
coronal iron abundances measured in the data and in the simula- 
tions may at first seem surprising, since it was shown in Paper I (see 
Fig. 12 of that study) that the hot CGM of the simulated galaxies 
is dominated by metal-poor gas accreted from the IGM. To under- 
stand how these findings are reconciled, we probe the enrichment 
of the simulated coronae in greater detail in this section. 

We remarked earlier, upon inspection of the surface density 
and surface brightness profiles presented in Fig. [2] that luminosity- 
weighted X-ray measurements provide an unrepresentative view of 
the state of hot CGM, because the X-ray emission does not trace 
the hot gas mass in a simple fashion. The effect of this bias on the 
inferred metallicity of the hot CGM is highlighted in Fig. H] where 
we plot the mass-weighted equivalent of Fig. [3] for the simulated 
galaxy sample. Each galaxy is represented by a filled circle whose 
colour encodes its halo mass, A/200, and a clear correlation between 
Lx and halo mass is evident (see also Fig. 6 of Paper I). Low lu- 
minosity galaxies exhibit a slightly greater median iron abundance 
and a larger scatter than their higher luminosity counterparts, be- 
cause the hot CGM associated with these systems is comprised of 
a greater fraction of gas recycled from evolved stars than is the 
case for more massive galaxies, as we discuss below. However, 
the median value of Zp e is typically < 0.1 Z across the range 
of luminosities probed here. This contrasts with the Lx-weighted 
measurement, which is near-solar across five decades in Lx- The 
colour coding reveals the strong correlation between Lx and halo 
mass (see also Fig. 6 of Paper I). For comparison, the inset panel 
shows the luminosity-weighed measurements from which the me- 
dian and scatter shown in Fig. [3] were obtained. 
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Figure 4. The mass-weighted Zp c — Lx relation of GIMIC galaxies. Galax- 
ies are represented by circles whose colour encodes the galaxy's halo mass, 
whilst the solid line, densely-hatched and sparsely-hatched regions show 
the median (hot gas) mass-weighted Z-p c , and its la and 2<r scatter, re- 
spectively. For reference, luminosity-weighted measurements (as in the left- 
hand panel of Fig. 15} are plotted in the inset. Although the Lx-weighted 
iron abundances are consistent with the near-solar values inferred from X- 
ray data, the hot CGM is metal-poor in a mass-weighted sense. Moreover, 
X-ray luminosity can be seen to correlate strongly with halo mass. These 
findings are expected if the hot CGM is established primarily via the accre- 
tion of gas from the IGM. 



The hot CGM of a typical ~L* galaxy in our simulation sam- 
ple at z — is therefore metal poor, indicating that the hot gas 
is not comprised primarily of the products of stellar evolution. It 
is possible to check this explicitly, since the gas mass 'recycled' 
by the stellar populations represented by star particles is shared 
amongst their neighboring SPH particles. Any difference between 
an SPH particle's mass at z = and the start of the simula- 
tion is therefore due solely to enrichment by stellar evolution. In 
Fig. B] we separate the simulated sample by halo mass into bins 
of Alog 10 M2oo[Mq] = 0.25. The left-hand panel shows the 
spherically-averaged radial profile of the median recycled mass 
fraction of particles comprising the hot CGM and, for reference, 
we plot in the inset the equivalent profiles for the entire CGM, 
irrespective of its thermal state. There is a strong radial depen- 
dence, such that the products of stellar evolution are, unsurpris- 
ingly, found close to the stars themselves, near the halo centre (met- 
als are, nonetheless, transported to large radius, as we discuss in 
§ |4.3[ l. The overall fraction of the hot CGM contributed by recy- 
cling is greatest for the lowest mass systems, as can also be seen 
in Fig. H] for two reasons. Firstly, the characteristic virial tempera- 
ture of these systems is low, so a smaller fraction of their accreted 
gas is X-ray luminous (and hence becomes part of the 'hot CGM'). 
Secondly, the baryon fraction of low mass haloes is lower than that 
of haloes hosting brighter galaxies because gas is more efficiently 
ejected by winds in the formeij The metal mass injected into the 



5 Whilst all of the galaxies of the HB06 and Ath07 samples are included 
in the 2MASS Extended Source Catalogue, the integrated fluxes of these 
nearby galaxies are recorded inaccurately, because of their proximity to 
scan edges. The 2MASS Large Galaxy Atlas accounts for this inaccuracy 
by combining adjoining scans. 



6 The baryon fraction of GIMIC galaxies is a strong function of halo po- 
tential. At z = 0, the typical baryon fraction increases smoothly from 10 
percent to 90 percent of the cosmic fraction for galaxies with circular veloc- 
ities between 100 < v c < 400 km s _1 (see Fig. 12 of C09), and exhibits 
significant scatter. This value is particularly sensitive to the details of the 
adopted feedback prescription (e.g. Haas et al.|2012), and it is likely that 
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Figure 5. Left: Spherically-averaged radial profiles of the fraction, by mass, of the hot CGM contributed by gas recycled by stellar evolution. Galaxies have 
been binned by Af200 (see legend). The inset shows the contribution to all halo gas (ISM + CGM), irrespective of thermal phase (note the different y-axis 
ranges). The contribution of recycled gas to the hot CGM is less than 10 percent at all radii (though the contribution to the ISM and cold CGM phases can 
be several times this). The small fraction of hot recycled gas is largely confined to small radii, close to the stars that produce it. Right: As left, but showing 
the fraction of soft X-ray luminosity produced by recycled gas. In spite of recycled gas contributing a small fraction of the hot CGM mass, it dominates the 
diffuse soft X-ray emission of the galaxy at all radii. This leads to the marked difference between the mass-weighted and luminosity-weighted metallicity seen 
in Fig.H] Note the markedly different scales on the y-axes of the three plots. 



hot CGM by stellar evolution is therefore diluted by a smaller mass 
of metal-poor gas accreted directly from the IGM. 

The overall contribution of recycled gas to the hot CGM is, 
nonetheless, small at all radii for haloes of any mass: even near 
the halo centre, recycled gas comprises less than 10 percent of the 
hot CGM (as also shown in Paper I), and integrated over the en- 
tire hot CGM the contribution is typically only a few percent. By 
contrast, as shown by the inset plot, recycled gas contributes a sig- 
nificant (but not dominant) fraction when considering all gas bound 
to the halo (i.e. ISM + all CGM gas, irrespective of thermal phase). 
Therefore a significant fraction of metals synthesised by the stars 
comprising today's i, galaxies have either been ejected into the 
IGM (see also |Ciotti et al.|1991} , or recycled into a cool phase (i.e. 
the ISM or cold CGM) from which they can become locked up in 
subsequent generations of stars (e.g. Leitner & Kravtsov 201 1 ). 

In the right-hand panel of Fig. B] we show the fraction of the 
diffuse soft X-ray luminosity produced by recycled gas. We com- 
pute the contribution from recycled gas on a per-particle basis, by 
equating the ratio Lx,rec/£x,tot to that of the specific soft X-ray 
emissivities, A rcc /Atot = 1 — A pr i/A to t, where A pr i is the specific 
emissivity of the gas packet assuming it has primordial element 
abundances. As discussed by Wiersma et al. (2009b), this approach 
is accurate as long as the contribution of heavy elements to the free 
electron density is a small fraction of the totalj The emissivities 
are computed as per Eq.|2] summing over 1 1 elements for each par- 
ticle, albeit with the element abundances (Zk in Eq.[2j of elements 
heavier than Helium fixed at zero in the case of A pr i. 

In contrast to the mass fraction shown in the left-hand panel, 
the luminosity fraction due to recycled gas dominates at all radii. 
This is because the specific emissivity of a ~ O.lkeV plasma 



the baryon fraction of massive galaxies here is overestimated owing to the 
omission of efficient feedback in massive haloes, such as that from AGN. 
7 For gas of solar composition, approximately one percent of electrons are 
contributed by metal ions. 



is extremely sensitive to the presence of metal ions. The effect 
is most dramatic in the case of low mass galaxies, which retain 
a smaller fraction of their baryons in the form of an accreted, 
quasi-hydrostatic corona (see C09), and in which galactic winds 
are more efficient at transporting enriched gas from the galaxy into 
the CGM and IGM in the GIMIC simulations. In addition, the bias 
between the recycled mass and luminosity fractions is expected to 
be largest for low mass galaxies, because the offset between the 
halo virial temperature, and the temperature at which the emissiv- 
ity of collisionally-excited iron peaks, is greatest for these systems 
l |van de Voort & Schaye|20T2) . Nonetheless, the emissivity of the 
hot CGM associated with even the most massive galaxies in our 
sample is dominated by recycled gas. Since the hot CGM of typ- 
ical galaxies remains extremely challenging to detect in emission, 
we infer from our simulations that the presence of metals in the hot 
CGM is a crucial aspect of its 'observability', since they exhibit a 
far higher specific emissivity than hydrogen and helium ions. 

It is important to clarify at this juncture what is meant by 're- 
cycled gas'. In the context of Fig.B] it would be more accurate to 
use the term 'recycled ions': the mass of the plasma is dominated 
by ions (as opposed to electrons), and when computing the plasma 
luminosity, we sum over individual ionic species. Physically, how- 
ever, the line emission is the result of collisions between ions and 
electrons (or, in the case of the Bremsstrahlung radiation that dom- 
inates more energetic plasmas, the result of the deflection of elec- 
trons by ions). Whilst the metal ions that dominate the emissivity of 
the plasma result from internal, stellar evolution processes, the left- 
hand panel of Fig. pj indicates that the mass density, and therefore 
the electron density, of the hot CGM is dominated by the metal- 
poor plasma accreted from the IGM. We therefore conclude that 
the observability (in emission) of the hot CGM of typical galax- 
ies requires both the outward transport of heavy elements from the 
central galaxy, and the accretion of gas from the IGM to compress 
and heat the recycled gas, and to provide electrons that can excite 
the metal ions. 
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Figure 6. Spherically-averaged radial profile of the mass-weighted (solid 
lines) and Lx-weighted (dotted lines) iron mass fraction of the hot CGM of 
GIMIC galaxies, binned by M200 as m Fig-P] The profiles of haloes of all 
masses exhibit a negative gradient. Mildly super-solar abundances are ex- 
hibited at the halo centre; the mass-weighted profiles drop to a few percent 
of solar at the virial radius, whilst the Lx -weighted profiles tend to a tenth 
of solar, highlighting the presence of a bias in the metallicity measure at 
large radii. The inset panel shows mean mass-weighted iron mass fraction 
obtained by integrating over all hot CGM gas within a radius r. See text for 
a discussion of these biases. 



3.3 Iron abundance mass profiles 

The declining contribution of gas recycled by stellar evolution 
to the mass of hot CGM at large radii is reflected in the radial 
profile of its mass-weighted (Eq.[3j iron abundance, Zf b , plotted 
with solid lines in Fig. [6] The strong abundance gradients reflect 
the centrally-concentrated Zf c distribution visible in the top-right 
panel of Fig.fT] The central iron abundance is typically super-solar, 
but drops to a few percent of Zf c ,q at the virial radius (r2on)- Such a 
strong metallicity gradient makes estimates of the hot CGM metal- 
licity dependent on the radius within which it is measured. To illus- 
trate this, we show in the inset plot the mass-weighted iron abun- 
dance recovered when integrating over the gas enclosed within a 
given radius, Zp c (< r). Within 0.1 r2oo (approximately r c for L* 
galaxies), the typical mass-weighted iron abundance is mildly sub- 
solar, but drops to ~ 0.1 Zq if all hot CGM gas within r2oo is 
considered. 

Our discussion so far has focussed on the fact that X-ray 
measurements of iron abundance are higher than the 'true' mass- 
weighted measurement because the strongly centrally-peaked sur- 
face brightness profile of the hot CGM strongly biases the mea- 
surement towards the bright, metal-rich centre of the X-ray corona. 
However, we might reasonably also expect a bias between mass- 
and Lx-weighted iron abundance measurements at fixed radius, 
because of i) the sensitivity of the plasma emissivity to the pres- 
ence of metal ions, and ii) the tendency of metals to clump (see e.g. 



Simionescu et al.|201 Tj Nagai & Lau|201 if | van de Voort & Schaye 



2012), either because they induce efficient cooling, and/or because 
they are associated with dense substructures. Clumping is indeed 
apparent in the stacked map of hot CGM iron abundance in Fig.fT] 
To assess the scale of this bias, we plot the Lx-weighted radial pro- 
file of Zp c in Fig.lSlwith dotted lines. At small radii (r < 0.2r2oo), 



the mass- and Lx-weighted measurements are similar, but at larger 
radii the mass-weighted iron abundance exhibits a steeper gradi- 
ent than the Lx-weighted measurement, such that at r2oo, the lat- 
ter is a factor of ~ 3 greater than the former. We conclude then 
that the effect is significant, but insufficient to explain the overall 
offset between mass-weighted and Lx-weighted measurements of 
Zp,, when integrated over the entire hot CGM. The combination 
of a centrally-peaked X-ray surface brightness profile and a strong 
negative gradient in radial metallicity in the hot CGM is stronger 
than the bias at fixed radius. 

In summary, it is clear that the 'true' (i.e. mass-weighted) iron 
abundance of the hot CGM is significantly lower than the value 
recovered from an Lx-weighted measurement. The marked dif- 
ference between the mass-weighted and Lx-weighted iron abun- 
dances reconciles the inference of near-solar iron abundances from 
X-ray spectroscopy of ~L* galaxies with models of galaxy for- 
mation in which the hot CGM is predominantly accreted from the 
IGM. The unavoidable Lx-weighting of X-ray data biases the mea- 
surements towards the hot, dense and metal-rich gas (with the high- 
est surface brightness) lying close to the galaxy. A mild bias be- 
tween mass- and Lx-weighted measures exists at fixed radius be- 
cause of the sensitivity of specific X-ray emissivity to the presence 
of metal ions. By mass, most hot CGM gas resides at large radius, 
where there is a greater mass of metal-poor gas, accreted from the 
IGM, to dilute transported metals. The rapid decline of the X-ray 
surface brightness of the hot CGM with galactocentric radius seen 
in Fig.fTJhighlights why radial metallicity profiles are so challeng- 
ing to determine from X-ray emission line spectroscopy. We spec- 
ulate that the best prospect for constraining mass-weighted abun- 
dance profiles with existing X-ray instrumentation is to use stacked 
emission maps or absorption spectroscopy (e.g. pang et al.| 2010), 
and to compare these to synthetic observations created with hydro- 
dynamical simulations. 



4 DISSECTING THE ENRICHMENT OF GALACTIC 
CORONAE 

Having established in §15] that the iron abundance of the hot CGM 
associated with >L* galaxies in the GIMIC simulations is compat- 
ible with observational constraints, we turn to a detailed analysis 
of how the enrichment proceeds over cosmic time. We start by ex- 
amining the location of enriched gas particles at the epoch of en- 
richment, and ask whether they were bound to the same central 
galaxy they are bound to at 2 = 0, a satellite galaxy of the central 
galaxy, another galaxy altogether, or bound to no galaxy at all (i.e. 
the IGM). We also explore the radial migration of hot, enriched gas 
throughout the assembly history of the galaxies. We subsequently 
investigate whether the enrichment of the hot CGM was dominated 
by prompt Type II SNe detonations or by the slow release of metals 
(on 10 9 -year timescales) by Type la SNe, AGB stars, etc. 



4.1 Is the hot CGM enriched internally or externally? 

To assess the relative significance of the reservoirs that contribute 
metals to the hot CGM of L* galaxies at Z — 0, we adopt a similar 
methodology to Font et al. ( 201 1 ), who explored the contribution of 
various processes to the formation and assembly of stellar haloes. 
For each galaxy in our GIMIC sample, we separate the enrichment 
history of its coronal gas into four categories: 

• in-situ: an enrichment event is flagged as 'in-situ' if, at the 
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time of enrichment, the recipient gas particle was bound to the most 
massive subhalo of the most massive progenitor (MMP) FoF group 
of the galaxy at z = 0. 

• satellite of MMP: the event is flagged as 'satellite enrichment' 
if, at the time of enrichment, the recipient gas particles was bound 
to any other subhalo of the MMP. 

• external galaxy: the event is flagged as 'external enrichment' 
if, at the time of enrichment, the recipient gas particle was bound 
to any subhalo of any FoF group other than the MMP. 

• smooth: the event is flagged as 'smooth' if, at the time of en- 
richment, the recipient gas particle was not part of any FoF group. 

The categories are unambiguous insofar that they account for all 
enrichment events and are mutually exclusive. 

The procedure by which we identify the MMP and classify en- 
richment events is as follows. For a given central galaxy at z — 0, 
we select all dark matter particles within r200 belonging to the dom- 
inant (most massive) subhalo. We use the unique IDs assigned to 
these particles to identify the FoF group in previous snapshots that 
contains the greatest number of these particles. This FoF group is 
said to be the MMP and the most massive subhalo of that FoF group 
is assumed to be the MMP of the dominant subhalo of the system 
at z — 0, Each of the gas particles that make up the hot CGM 
at z — is tracked back in time through all previous simulation 
outputs ('snapshots'), and the iron mass added to the particle is as- 
signed to one of the four categories above, based on which, if any, 
FoF group and subhalo the gas particle belonged to at that time. 

Each GIMIC volume has approximately 60 snapshots, span- 
ning the redshift range z — 20 to z — 0. The time of each enrich- 
ment event will not, in general, correspond to the redshift for which 
the snapshot data were written. We therefore must use the snap- 
shot that is closest in time to identify where the gas particle was 
enriched. This could, in certain circumstances, lead to a misiden- 
tification of the location of the enrichment (i.e. to which, if any, 
subhalo it belonged when it was enriched). However, we have ex- 
plicitly checked that this is not an important issue by re-running our 
tracing algorithm using only every second snapshot, and find that 
the in-situ enrichment fraction recovered in this case does not dif- 
fer significantly; we therefore conclude that the redshift sampling 
of our snapshot data is sufficient to deduce whether a given par- 
ticle was enriched 'internally' (i.e. in-situ) or 'externally' (i.e. in 
a satellite, an external galaxy, or the IGM). We trace each of the 
z — coronal gas particles back in time as far as z = 10 to deter- 
mine their enrichment category. We ignore the very small fraction 
of metal mass synthesised at z > 10. 

In Fig. [7] we plot the fractional contribution of these categories 
to the overall iron mass fraction radial profile, for the ~ 350 galax- 
ies with halo mass in the range 10 1200 < A# 2 oo/M© < 10 12,50 . 
Within the virial radius in-situ enrichment dominates, indicating 
that the majority of coronal metals are synthesised by, and trans- 
ported from, stars comprising the central galaxy. At large radii (~ 
?"20o), the 'external galaxy' category becomes significant (but not 
dominant), indicating that coronal outskirts are, in part, enriched by 
the transfer of gas from external galaxies to the MMP, for example 
via ram-pressure stripping or entrainment within supernova-driven 
(SN-driven) outflows. It is clear that the contribution of external 
galaxies to the metal mass of hot gas confined by the potential of 
haloes will be a function of halo mass, since the central galaxy of 
galaxy groups and clusters hosts a decreasing fraction of the overall 
stellar mass in more massive haloes (Lin & Mohr 2004). However, 
even in haloes as massive as ~ 10 1 M© the central galaxy typi- 
cally comprises over 50 percent of galactic stellar mass bound to 
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Figure 7. Spherically-averaged mass-weighted radial profiles of the coronal 
iron mass fraction contributed by the categories of enrichment described in 
§ |4.1| in-situ (blue), within an external galaxy (red), within a satellite of 
the central galaxy (green), and external to any FoF halo ('smooth', cyan). 
Within T2oo, in-situ enrichment dominates, indicating that the majority of 
coronal metals are transported from the stars comprising the central galaxy. 
At large radii, where relatively few in-situ metals reside, metals aqequired 
from external galaxies (e.g. by winds or ram-pressure stripping) become 
significant, but not dominant. Satellites of the main progenitor of the galaxy 
contribute only a small fraction of the total metal content of the hot CGM. 



the halo, so the curves presented in Fig.lTlare unlikely to be partic- 
ularly sensitive to M2oo- 

The enrichment of gas associated with satellites of the MMP 
contributes, on average, only a small fraction of the total iron mass 
of the simulated coronae. This is, in part, due to our sample selec- 
tion criterion that the main subhalo of a FoF group should account 
for at least 90 percent of the virial mass, a criterion we imposed 
to ensure that we study the hot coronal gas associated with galax- 
ies rather than major mergers or galaxy groups. However, we have 
confirmed that this bias is small by relaxing this criterion to require 
that M su b > 0.5Mf o f, which produces very similar results. 

The contribution to the overall fraction of coronal iron from 
intergalactic gas ('smooth') can be viewed as an uncertainty esti- 
mate on the plotted fractions. Since GIMIC assigns the metals syn- 
thesised by star particles to their neighbouring gas particles (i.e. 
those within an SPH smoothing kernel), in the limit of infinite 
resolution nearly all gas particles should be enriched as part of a 
FoF (b — 0.2) group (the fraction of truly intergalactic stars is 
very small). However, the finite particle resolution of the simulation 
means that it is possible for the smoothing kernel of star particles 
to enclose a small fraction of unbound particles. Fortunately, this 
issue is most problematic in ill-resolved, low-mass galaxies, which 
synthesise only a small fraction of the cosmic heavy elements, as 
confirmed by Fig. [7] We therefore neglect this component in sub- 
sequent plots. We note that our findings are qualitatively consistent 
with those of Shen et al. (2012), who simulated the formation of a 
massive galaxy to z — 3 and found that the central galaxy was the 
dominant source of circumgalactic metals, with only ~ 10 percent 
contributed by satellite companions at that epoch. 

In addition to ascertaining the mode by which gas is enriched, 
we can obtain a simple quantitative diagnostic that elucidates where 
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Figure 8. Spherically-averaged mass-weighted radial profile of the iron en- 
richment weighted radius (Eq. [31 of coronal gas at z = 0, split by the 
enrichment categories described in § |4.1| The diagonal dotted line marks 
the locus r-pc = r(z = 0), denoting no radial migration. The mode of iron 
synthesis correlates strongly with the galactocentric radius at which it was 
synthesised, and both inwards and outwards metal transport are significant. 
Overall, outward transport dominates, since most coronal iron is produced 
in-situ and re-distributed by galaxy-scale winds and convection. 




Figure 9. Spherically-averaged mass-weighted radial profile of the iron en- 
richment weighted redshift (Eq. 17) of coronal gas at z = 0, split by the 
enrichment categories described in § |4.1| The enrichment categories differ 
markedly, with enrichment by external galaxies occuring at relatively early 
times and in-situ enrichment being dominant at late epochs. All categories 
exhibit a trend such that the metals located farthest from the galactic cen- 
tre were enriched at the earlist epochs. We interpret this as a signature of 
a travel-time, such that the metals synthesised earliest have been in transit 
from the galactic centre for the longest time. 



coronal gas was enriched by defining a quantity for each gas parti- 
cle that we term the 'iron enrichment weighted radius', rp e : 



rp c = 



Y^ z Ami, Fe ri 

Y^i Am j>Fe 



(6) 



where i sums over all enrichment events, Ami,F e is the growth in 
iron mass of the gas particle per event, and n is the radial distance 
(in physical coordinates) of the enriched gas particle from the halo 
centre at the time of enrichment (see also Wiersma et al. 2010} |Shen"l 
|et al.|20 l2). This quantity is not followed explicitly on a timestep- 
by-timestep basis in GIMIC, hence we compute it in post-processing 
from snapshots. 

In Fig. [8] we plot ?~fc as a function of z — radius, binned by 
the enrichment categories defined above. Iron synthesised within 
external galaxies (red) was, by definition, synthesised beyond the 
virial radius of the MMP, but later became bound to the potential of 
the central galaxy-halo system. This can occur as a result of metal 
transport due to winds driven by the external galaxy (e.g. |Shen| 
|et al.| |2012), or, if the external galaxy subsequently merges with 
the MMP, due to ram pressure stripping of the external galaxy's 
CGM as it interacts with the hot CGM of the MMP ( [Abadi et al] 
1999; Bahe et al. 2012). Externally synthesised iron exhibits a mild 
correlation between the tf c and r(z = 0), such that the iron syn- 
thesised closer to the coronal centre resides closer to the galactic 
centre at z — 0. 

The dominant in-situ (blue) mode of iron synthesis exhibits 
a significant correlation between rp c and r(z — 0). Interestingly, 
this curve crosses the dotted line in Fig. [8] that marks the locus 
rF C = r(z = 0), and hence the in-situ iron is characterised by 
both inward and outward transport. The inward transport is con- 
fined to iron synthesised within ~ 0.1 of the present day virial ra- 
dius of galaxies, which roughly corresponds to the scale length of 
star forming galaxy discs: this iron was likely therefore synthesised 



in the outskirts of galaxy discs and subsequently migrated radially 
in the disc plane. However, most in-situ iron is transported away 
from the central galaxy during its assembly history. We therefore 
conclude that whilst both internal and external proceses contribute 
to the enrichment of the hot CGM associated with L* galaxies at 
z = 0, it is the internal, in-situ category that dominates. 



4.2 When were coronae enriched? 

In analogy with the iron mass weighted radius defined in Eqn.[6] we 
can we can obtain a simple quantitative diagnostic that elucidates 
•when coronal gas was enriched by defining an 'iron mass weighted 
redshift', zp c , for each particle (see also Wiersma et al. 20101: 

J2i Ami,Fe2> 



ZFb = 



£\ Am,,Fc 



(7) 



where, i once again represents an enrichment event, Ami.p c is the 
growth in iron mass of the gas particle per event, and Zi is the 
epoch of enrichment. As for rp a , we compute this quantity in post 
processing and, in analogy with Fig. [8] plot zp c as a function of 
z = radius, binned by the enrichment categories defined in § |4.1| 
in Fig. 15] The enrichment categories also exhibit significant differ- 
ences in this diagnostic. 

Iron contributed by external galaxies was synthesised at early 
times (typically 1 < z < 3), and a mild correlation between 
r(z — 0) and zp c is apparent, such that the externally-synthesised 
iron close to the coronal centre at z — was synthesised more re- 
cently than that in the coronal outskirts. The in-situ iron exhibits a 
strong correlation between riz = 0) and zp c , such that the most- 
recently synthesised iron is closest to halo centre. In-situ iron found 
within 0.1r2oo was typically synthesised within the last 3 — 4 Gyr, 
whilst the in-situ iron at ~ r2oo was synthesised 7—8 Gyr ago. This 
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trend is a clear signature of the operation of metal transport, insofar 
that the iron synthesised at the earliest epochs has had the longest 
period of time to be transported to large radii. Interestingly, qualita- 
tively similar trends operate for the satellite and external galaxy cat- 
egories, indicating that the majority of metals synthesised in satel- 
lite and external galaxies are likely deposited close to the halo cen- 
tre, and transported outwards in the same fashion as the in-situ iron. 
This result is in qualitative agreement with Wiersma et al.| ( [20T0) , 
who found that metals residing in low density gas were typically 
ejected at early times, by low mass haloes. 



4.3 How are the metals transported? 

Having established that the majority of metals in the hot CGM of 
today's L* galaxies were synthesised at the halo centre and trans- 
ported outwards during the assembly history of the galaxy, it is 
interesting to consider briefly how the transport proceeds. The en- 
trainment of metals in galaxy winds is a natural mechanism to con- 
sider, since the winds originate at the galactic centre and propagate 
into the CGM. Moreover, metals entrained in outflows are readily 
observed at both high and low redshift (e.g. Heckman et al. 1990, 
|Pettini et al.|2001||VeTileux et al.|2005| >. The importance of metal 
transport by winds as a means to enrich circum- and intergalactic 
gas was recently studied in detail by Wiersma et al. (2011), who 
used cosmological hydrodynamical simulations to show that winds 
are the primary mechanism by which the cold-warm IGM is en- 
riched. 

However, we demonstrated in Paper I that the hot CGM of 
L+ galaxies in the GIMIC simulations is dominated at late-times 
by quasi-hydrostatic gas that has never been part of the ISM (and 
hence has never been ejected in a galactic fountain). It is there- 
fore unlikely that outflows are the dominant metal transport mech- 
anism at late times (z < 1) in GIMIC. There is, however, another 
mechanism by which the metals can be transported: convection. 
|McCarthy et al.| ( |2007| l demonstrate that hydrostatic haloes grow 
in an inside-out fashion, such that the orbital kinetic energy asso- 
ciated with infalling gas is thermalised within the hot CGM of the 
main progenitor. The entropy of the infalling gas itself is largely 
unaffected, so it sinks to the halo centre and thus steadily convects 
the higher entropy and enriched gas, initially at the halo centre, to 
larger radii. It follows naturally that the metals synthesised earliest 
will reside furthest from the halo centre. Although some fraction of 
the hot CGM will flow inwards due to radiative cooling, a signifi- 
cant fraction of the gas has a long cooling time and will therefore 
convect as the halo is assembled. The dominance of this mecha- 
nism on timescale grounds appears plausible; transporting metals 
a distance of 250 kpc (the typical r2oo at z — for a galaxy in 
our simulation sample) in a period of 8 Gyr requires a net outflow 
velocity of only ~ 30 km s _1 , roughly an order of magnitiude 
lower than the characteristic velocity of SN-driven outflows as they 
emerge from the ISM. 



4.4 Which mechanism dominates iron synthesis? 

Having analysed the radial distribution of iron in the hot CGM, and 
the nature of its release into, and transport within, the hot CGM, 
we now turn to the source of the iron. We consider the synthesis by 
Type II SNe that occurs promptly (< 30 Myr) after star formation 
episodes, and the synthesis by Type la SNe that is delayed with re- 
spect to star formation. As described by Wiersma et al. (2009bl, the 
chemodynamics implementation adopted by GIMIC also accounts 
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Figure 10. Spherically-averaged radial profile of the mass fraction of iron 
within the hot CGM, synthesised by Type la SNe. The median iron mass 
fraction synthesised by Type la SNe within 7*200 for each halo mass bin 
is quoted next to the legend. Iron associated with the hot CGM of GIMIC 
galaxies is predominantly synthesised by Type la SNe; this is in spite 
of our adopted yields and Type la SN rates being such that an evolved 
(tagc > 10 9 Gyr) SSP with a Chabrier IMF tends to a near-equipartition 
of iron production from Type II SNe and Type la SNe. We infer that iron 
synthesised by Type II SNe is preferentially i) entrained in outflows and 
ejected from the system, or ii) locked up in secondary generations of stars. 



for the fact that, whilst AGB stars do not synthesise iron, they do 
'lock-up' significant quantities of iron when an SSP forms from 
pre-enriched gas. This iron is then gradually released by AGB en- 
velope shedding on ~ Gyr timescales. 

In Fig. [TO] we split the profiles presented in Fig. [6] into con- 
tributions from the two types of SNe, and plot the spherically- 
averaged radial profile of the Type Ia-to-total iron mass fraction. 
The radially-integrated (0 < r < r-2oo) iron mass fraction synthe- 
sised by Type la SNe is quoted for each halo mass bin in the legend. 
Both types of SN contribute at all radii for galaxies of all masses. 
The approximate equipartition of iron synthesised by each SN type 
follows from the behaviour of the stellar populations comprising 
the galaxy associated with each corona; adopting the yields and 
Type la SNe rates presented by Wiersma et al. (2009b), an evolved 
(iagc > 10 9 yr) SSP with a Chabrier IMF tends to a broad equipar- 
tition of the iron mass fraction produced by each type of SN, par- 
ticularly if the SSP has sub-solar metallicity. However, we caution 
that this statement is sensitive to the assumed cosmic Type la rata*] 
and the iron yield of Type la and II SNe. 

There is, however, a mild but significant preference for iron 
synthesised by Type la SNe over that from Type II SNe in all but 
the least massive galaxies. We deduce two reasons why iron syn- 
thesised by Type II SNe is inhibited from residing in the hot CGM 
of galaxies at late cosmic epochs: firstly, the iron sythesised by 
Type II SNe is preferentially entrained in the energetic outflows 
that these events trigger, enabling a significant fraction of iron syn- 



8 The SN la delay function is ill-constrained by theory and observation. 
The GIMIC simulations adopt an e-folding delay function with a character- 
istic timescale of r = 3.3 Gyr. 
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thesised in this fashion to be ejected into the IGM at early epochsj 
Secondly, because Type II SNe synthesise iron promptly, any iron 
from Type II SNe that is not ejected from the galaxy's potential 
is deposited into the dense ISM, where the enrichment products 
trigger efficient radiative cooling and readily become locked-up in 
secondary generations of stars. This is reflected in the abundance 
patterns of stars in the simulated galaxies: the fraction of Type II 
SN-synthesised iron in the stellar component of galaxies compris- 
ing our sample is typically 70-80 percent (see also |Font et al.|20lT| >. 
However we caution that this is likely a mild overestimate of the 
true fraction, since the inefficiency of feedback in the most massive 
galaxies comprising our sample will lead to too high a fraction of 
Type II SN-synthesised iron being unable to escape the confines of 
the ISM. 

The trend with halo mass, such that the hot CGM of the most 
massive haloes is more strongly dominated by iron synthesised by 
Type la SNe (^F C ,SNia/^Fc,tot ~ 60 percent) can be interpreted in 
terms of the different star formation histories of the galaxies associ- 
ated with these coronae. As shown in C09 (see Fig. 6 therein), more 
massive galaxies in GIMIC complete a greater fraction of their star 
formation at high-redshift, enabling a greater fraction of their stel- 
lar populations to evolve to the stage where a significant fraction of 
their ejected iron is contributed by Type la SNe. 

The coronal iron mass fraction from Type la SNe in GIMIC is 
slightly lower than the 70 — 90 percent estimated by HB06 for their 
sample of ellipticals, whose masses and X-ray luminosities are con- 
sistent with the most massive galaxies in the GIMIC sample. How- 
ever, comparison of this result is not entirely instructive; we have 
not limited the GIMIC sample to 'elliptical' galaxies, and the lim- 
ited X-ray spectroscopy available for disc galaxies with high star 
formation rates unsurprisingly indicate abundance patterns more 
consistent with Type II SN yields (Riching s et al.|2010[ l.[Tang & 
Wang (2010) argue that the observable iron content of the hot CGM 
might in fact provide an unpresentative view of the total mass of 
iron synthesised by Type la SNe, because the iron is likely to be 
entrained in very hot, tenuous ejecta with a low specific soft X-ray 
emissivity and a large buoyancy that results in its rapid transport to 
large radii. Perhaps more importantly, however, SN nucleosynthe- 
sis yields remain uncertain at the factor ~ 2 level, and the cosmic 
Type la SN rate remains poorly constrained by both theory and ob- 
servation. The simulations could remain entirely compatible with 
the cosmic Type la SN rate using double the adopted rate, and this 
would unquestionably improve the correspondence between the ob- 
served and predicted Type la iron fraction in the hot CGM. For this 
reason, we caution that comparison of the simulations with the cur- 
rently available observational constraints is not particularly instruc- 
tive. 



5 SUMMARY AND DISCUSSION 

We have investigated the enrichment of the hot circumgalactic 
medium (CGM) associated with ~L* (10 10 < M* < 10 115 M ) 
galaxies at z = 0, using the GALAXIES-lNTERGALACTIC 
Medium Interaction Calculation (gimic; |Crain et al7| 
2009). These simulations produce a large sample (n = 617) of 
well-resolved galaxies at z — 0, with a distribution of morpholo- 
gies that broadly corresponds to that observed in the local Universe. 



9 This is clearly a significant process, since the mass of gas ejected by 
a galaxy's stars over cosmic history is typically comparable to, or greater 
than, the mass of gas associated with the galaxy's potential at z = 0. 



We showed in |Crain et al.| {2010 Paper I) that the simulations re- 
produce the observed z = scalings of X-ray luminosity with K- 
band luminosity, star formation rate, and disc rotation velocity. The 
simulations also reproduce the observed z = stellar mass - rota- 
tion speed (or 'Tully-Fisher') and stellar mass - halo mass relations 
for 10 9 < M+/M Q < lO 10 ' 5 , but we caution t hat they still suf- 
fer from some overcooling for M* > 10 11 Mq i McCarthy et al 
2012b). The results from our simulations are summarised as fol- 
lows: 

• The hot CGM is dominated in a mass-weighted sense by 
metal-poor gas accreted from the intergalactic medium. Therefore 
the iron abundance of the hot CGM, when integrating over all gas 
within the virial radius, is typically < 0.1Zp e> Q. 

• Conversely, the X-ray emissivity of the hot CGM is dominated 
by line emission from collisionally excited metal ions that were 
synthesised in stars and transported out of the ISM. However, the 
excitation is caused by collisions with electrons; since heavy ele- 
ments make only a small contribution to the free electron budget, 
these electrons are primarily sourced from the metal-poor gas ac- 
creted from the IGM. 

• The iron found in the hot CGM of galaxies in the GIMIC simu- 
lations at z = was predominantly injected into gas that was, at the 
epoch of enrichment, already bound to the most massive progenitor 
of the galaxy, a process we term 'in-situ' enrichment. 

• We find a clear correlation between the galactocentric radius 
of metals at z = 0, and the redshift at which they were synthesised, 
such that those synthesised at earlier cosmic epochs reside farther 
from the galaxy. Enrichment therefore proceeds in an 'inside-out' 
fashion, requiring the transport of metals from the central galaxy 
into the CGM. We speculate that this transport is dominated by 
SNe-driven winds at early times, and by convection associated with 
the inside-out growth of the halo at later times, when the hot CGM 
is quasi-hydrostatic. 

• Inside-out enrichment establishes a strong negative iron abun- 
dance gradient in the hot CGM. The spherically-averaged radial 
profile of the iron abundance is typically super-solar at small radii 
(< 0.05r2oo)> an d rapidly declines to only a few percent of the so- 
lar iron abundance at 7-200, or to about ten percent of solar if the 
metallicity is weighted by the luminosity. . 

• The strong negative metallicity gradient established by inside- 
out enrichment concentrates the majority of the metal ions associ- 
ated with the hot CGM inside a relatively small radius. In contrast, 
following simple volumetric arguments, the majority of the mass 
of the hot CGM resides at large radius. Therefore, X-rays do not 
trace the distribution of the hot CGM in a simple fashion, rendering 
metallicity measurements of the hot CGM inferred from spatially 
unresolved X-ray spectroscopy biased towards the metal-rich gas 
close to the galaxy. 

• In addition, the sensitivity of plasma emissivity to the pres- 
ence of metal ions induces a bias such that the luminosity- weighted 
metallicity of the hot CGM at a fixed radius is typically elevated 
with respect to its mass-weighted metallicity at the same radius. 
This induces a bias between luminosity- and mass-weighted metal- 
licity measurements even when inferred from spatially-resolved X- 
ray spectroscopy. 

• Both biases lead to abundance measures that are elevated with 
respect to the 'true', mass-weighted, measure. Therefore, despite 
exhibiting a typical iron abundance of < 0.1Zp Oi Q at z = 0, 
the median Lx -weighted iron abundance of the hot CGM is ap- 
proximately Zf Ci 0, with a 2<r scatter of approximately one decade, 
across 5 decades in Lx- 
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The simulations yield a correlation between Zf c and Lx that 
is weaker than implied by the combination of the Hump hrey &| 
|Buote| < [2006[ > and |Athey| ((2007) datasets, but the recovered iron 
abundances are broadly consistent (within the 2a scatter) with these 
data. The precise nature and origin of this relation is clearly an 
avenue of interesting future theoretical and observational enquiry. 
Similarly, the median core radius of our hot CGM profiles is larger 
than the core sizes inferred from the small number of measurements 
inferred from deep exposures of individual galaxies, or stacked pro- 
files extracted from the ROSAT All-Sky Survey. As discussed in 
§ |2.7| we do not consider this a crucial shortcoming of our model, 
but note that the establishment of the surface density profile of cir- 
cumgalactic gas remains an unresolved astrophysical problem that 
is well suited for study with the next generation of hydrodynamical 
simulations. 

5.1 The interpretation of X-ray data in the context of galaxy 
formation theory 

As remarked in § |3.1[ a central aim of this study is to ascertain 
whether the general picture of the formation of the hot CGM that is 
predicted by galaxy formation theory, and seen in cosmological hy- 
drodynamical simulations, is compatible with constraints inferred 
from the most authoritative X-ray datasets presently available. 

The GIMIC simulations follow the hierarchical assembly of the 
galaxies and include detailed phenomenological models for the for- 
mation of stars, the effects of SN-driven feedback, and the chemical 
evolution of evolving stellar populations. Whilst the simulations 
are known to produce inaccurate results in certain regimes (e.g. 
C09, [McCarthy et al. 2012b) and necessarily offer an incomplete 
description of the formation of galaxies, they can be considered 
an authoritative testbed for confrontation with observational mea- 
surements, by virtue of the broad range of observable properties 
of local ~L* galaxies that they successfully reproduce (Fo nt et al.| 
[20TT]pcC~arthy et al.|2012a|b) . 

In Paper I and in this study, we establish that ~L* galaxies 
in GIMIC broadly reproduce the key observed X-ray scaling rela- 
tions of galaxies (i.e. with rotation velocity, stellar mass and star 
formation rate), and their inferred luminosity-weighted hot CGM 
iron abundance. Although the simulations are formally not fully 
consistent with the HB06 and Afh07 datasets, there is a broad cor- 
respondence between the two. This is a significant result, since the 
relatively high metallicity of the hot CGM of local galaxies has 
been perceived as a challenge to the prevailing model of galaxy 
formation, within which the hot CGM is assumed to form primarily 
via the accretion of metal-poor gas from the intergalactic medium 
throughout the assembly of the galaxy. 

Our simulations indicate that the hot CGM of L* galaxies 
is indeed dominated by gas accreted from the IGM. This picture 
appears similar to that advanced by simple analytic/semi-analytic 
models (e.g. |White & Rees|1978||White & Frenk|1991| >, but in fact 
differs in two important ways. Firstly, in contrast to the common 
assumption that the radial density profile of the hot CGM traces 
that of the underlying dark matter distribution, we showed in Paper 
I that SN-driven feedback lowers the normalisation of the density 
profile and acts to soften the central cusp into a core by preferen- 
tially ejecting low entropy gas from the halo centre. The reduction 
in central gas density lowers the luminosity of the halo by up to two 
orders of magnitude, because both line emission from collisional 
excitation and Bremsstrahlung radiation vary with the square of the 
gas density. 

Secondly, as shown here, the dynamically complex and time- 



dependent enrichment of the hot CGM by the central galaxy, its 
satellites and external galaxies, establishes a strong negative radial 
metallicity gradient. This contrasts with the assumption of com- 
plete, uniform mixing within the cooling radius adopted by WF91. 
Since the emissivity of the gas is particularly sensitive to the pres- 
ence of metal ions, this is an important difference that, like the 
density profile, has a marked effect on the luminosity and surface 
brightness profile of the hot CGM. Therefore, it is particularly en- 
couraging that, in contrast to WF91, GIMIC broadly reproduces the 
luminosity and surface brightness profiles of galaxies in the local 
Universe. 

This notwithstanding, the interpretation of diffuse X-ray emis- 
sion associated with local ~L* galaxies remains controversial: the 
most common interpretation in the literature is that the emission as- 
sociated with disc galaxies is a signature of hot plasma entrained in 
outflows driven by Type II SNe (e.g. Strickland et al. 2004; [Wang] 
|2005| |Tullmann et~aT1|2006| |Rasmussen et al.||2009| >, whilst that 
associated with ellipticals is due to the heating of cold ISM gas 
by its interaction with the ejecta of intermediate age (AGB) stars 
and Type la SNe (e.g. Mathews 1990; Ciot ti et al.|1991||Mathews| 
|& Brighenti|2003||Parr iott & Bregman 2008}. Whilst simple mod- 
els based on these premises can lead to suggestive links between 
the optical and X-ray properties of galaxies, such models are, by 
construction, incomplete, because they start from idealised initial 
conditions, and are often unphysical, for example when assuming 
that galaxies are hydrodynamically decoupled from the intergalac- 
tic medium. Moreover, such models neglect, by construction, the 
role of the CGM as a reservoir of fuel for on-going star formation 
in galaxies. As we have argued in this study, a more complete model 
of galaxy formation is necessary to study the X-ray emission of the 
hot CGM, because it appears that the emission is the result of both 
accretion from the IGM and the transport of heavy elements from 
the ISM into the hot CGM. 

The outstanding challenge is therefore to accomodate the in- 
terpretation of X-ray observations of galaxies within the context of 
modern hierarchical galaxy formation theory. As demonstrated in 
Paper I and in this study, many such observations are in fact re- 
produced by realistic, self-consistent hydrodynamical simulations 
that adopt cosmological initial conditions. For example, the broad 
correlations between Lx and AL (in disc galaxies) and AL (in el- 
lipticals) that provide the fundamental motivation for the idealised 
models described above, arise naturally in the GIMIC simulations 
because Lx, M* and AL all correlate with M200 for ~L* galaxies. 
It seems not unreasonable to expect, therefore, that improved mod- 
elling in this context will enable future simulations to reproduce 
a wide range of X-ray properties and scaling relations accurately. 
Moreover, as remarked in Paper I, the fundamental nature of the hot 
CGM of typical galaxies is, in principle, directly accessible via X- 
ray line diagnostics. Future X-ray observatories with the ability to 
obtain high-sensitivity spectroscopy at high spatial and spectral res- 
olution therefore offer a direct means to identify the fraction of hot 
CGM entrained in outflows or residing in quasi-hydrostatic reser- 
voirs. Observations with these facilities may thus verify or falsify 
the key role in the formation and evolution of galaxies that modern 
theoretical galaxy formation models assign to hot gaseous coronae. 
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APPENDIX A: NUMERICAL CONVERGENCE 

We next consider the effects of numerical resolution. Focussing on 
the Zpe-Lx scaling relations, we follow the methodology adopted 
by C09, |Font efaTlpOTT) and |McCarthy et al.| ( [20T2a) l, and com- 
pare the Zf b — I/x relation derived from the intermediate-resolution 
simulations with results from the high-resolution realisation of the 
—2a GIMIC region. The high resolution realisations have a factor 
of eight better mass resolution than the intermediate-resolution re- 
alisations (see § |2.1fr . 

In Fig. |A1| we compare the mass-weighted (left) and Lx- 
weighted (right) metallicity-luminosity scaling of the 74 galaxies 
satisfying our selection criteria (see § |2.3[ > in the high-resolution 
—2cr region (black points), with the median and scatter derived 
from the sample of 617 galaxies in all five intermediate-resolution 
regions (red hatching). The correspondence between the high- 
and intermediate-resolution samples is good, indicating that the 
Zf c — Lx relations (both mass-weighted and Lx-weighted) are nu- 
merically convergerged in the intermediate resolution simulations. 



APPENDIX B: METALLICITY SMOOTHING 



In § |2.6| we discussed the lack of implicit scalar quantity diffusion 
in SPH, and the consequent suppression of metal mixing that un- 
avoidably results from the finite sampling of the fluid. Here, we 
suppress the effects of suppressed mixing by kernel-weighted aver- 
aging element abundances over a given particles's SPH neighbours 
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Figure Al. The mass-weighted (left) and luminosity-weighted (right) Zp c — Lx relations defined by the 74 galaxies satisfying our selection criteria in the 
high-resolition —2a GIMIC simulation, overplotted on the median and (1, 2)<r-scatter (hatched regions) defined by the sample of 617 galaxies drawn from all 
five intermediate-resolution simulations. The relations described by the intermediate- and high-resolution simulations are numerically converged. 



when computing cooling rates, radiative luminosities, and coronal 
metallicities. 

This smoothing proced ure induces only a sm all error in the re- 
sulting total metal mass (see Wiersma et al. 2009b| 10 ] but does have 



the potential to bias the recovered metallicity of hot circumgalac- 
tic particles, because smoothed metallicities are computed using 
all SPH neighbour particles, irrespective of their thermal state. The 
smoothed element abundances of hot particles close to the disc- 
corona interface are therefore computed from both cold (i.e. ISM) 
and hot (coronal) particles. The maximum possible effect of such 
a potential bias can be determined by recomputing the metallic- 
ity and X-ray luminosity of each galaxy's hot CGM, using particle 
abundances instead of the SPH-smoothed abundances. 

In Fig. |Bl| we compare the mass-weighted (left-hand column) 
and Lx-weighted (right-hand column) Zp c — Lx relations derived 
from the particle (top row) and smoothed (bottom row) abundances. 
As in previous figures, each dot represents a simulated galaxy, 
whilst the densely- and sparsely-hatched regions represent the la 
and 2a scatter, respectively. 

It is clear that the total X-ray luminosity is insensitive to metal 
mixing uncertainties, since no systematic shift in Lx is observed. 
Similarly, the recovered mass-weighted metallicity is little changed 
in well-resolved coronae, since the smoothing broadly conserves 
metal mass. In the lowest luminosity systems, which are less well 
resolved than more massive counterparts, the coronal metallicity 
derived from the particle abundances is lower by a factor of < 2, 
indicating that the degree of biasing caused by the smoothing of 
metals from cold ISM is small. It is, for example, smaller than the 
uncertainty in our adopted metal yields. 

The effect on the luminosity-weighted metallicity is slightly 
larger: metallicities derived from particle abundances are lower by 
a factor of < 3. In contrast to the mass-weighted quantities, this 
result is relatively insensitive to the coronal luminosity, because the 
most luminous particles are always those closest to the disc-corona 



interface. Therefore, the extent of the largest possible bias resulting 
from our use of smoothed abundances is only marginally greater 
than the intrinsic uncertainty in our adopted metal yields, and does 
not impact upon our conclusions. 



APPENDIX C: PROPERTIES INFERRED BY HB06 & 
ATH07 FROM Chandra SPECTROSCOPY 

As described in § |3.1| we use the HB06 sample as our fiducial sam- 
ple and supplement this with unique galaxies from the Ath07 sam- 
ple, since these extend to lower X-ray luminosity. It is interesting 
briefly to examine the correspondence between the values of Zf c 
and Lx recovered by both studies from the archival Chandra spec- 
troscopy for the 22 galaxies that were included in both samples, in 
order to assess the impact of systematic uncertainties. 

In Fig. |Cl| we plot the soft X-ray luminosity (left) and the iron 
abundance (right) of hot circumglactic gas, reported by Ath07, as 
a function of that reported by HB06. In each case, the 1 : 1 corre- 
spondence is shown as a dotted line. Ath07 report systematically 
lower luminosities, with the faintest and brightest galaxies differ- 
ing by over a decade in Lx, whilst intermediate luminosity galaxies 
(Lx~ 10 42 erg s _1 ) are offset by ~ 0.5dex. The inferred metal- 
licities show no such systematic offset, and broadly trace the 1 : 1 
locus, albeit with broad scatter. 

The differences between the measurements highlight the diffi- 
culty of reducing X-ray spectra, and fitting plasma models to these 
data to infer luminosities and metallicities. However, the uncertain- 
ties apparent here do not affect our conclusions: the inferred metal- 
licities are not systematically offset, whilst the application of a sys- 
tematic offset to the luminosties reported by either HB06 or Ath07 
would merely require the data to be translated along the x-axis in 
the left-hand panel of Fig. [5] Since the simulated galaxies exhibit a 
flat Zf -Lx relation, such a translation would leave the simulations 
entirely compatible with the data. 



u Strict conservation is not achieved because metals are distributed with a 
gather, rather than scatter, approach 
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Figure Bl. Comparison of the mass-weighted (left-hand column) and Lx-weighted (right-hand column) Zp c — Lx relations recovered when adopting particle- 
based (top row) and smoothed (bottom row) element abundances. The different abundance measures produce differences in the mass-weighted relation only 
for low-luminosity systems that are poorly resolved. The luminosity weighted relation exhibits differences for all systems probed here, but only at the factor 
< 3 level. Our conclusions are therefore insensitive to the use of the particle-based or smoothed element abundances when measuring coronal metallicities 
and computing particle X-ray luminosities. 
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Figure CI. The X-ray luminosity (left) and iron abundance (right) inferred by Ath07 as a function of that inferred by HB06, for the 22 galaxies present in both 
samples. The names of these galaxies are listed in the left-hand panel. Luminosity uncertainties are quoted only by Ath()7. As in Fig. [3] we plot galaxies from 
HB06 with no formal upper limit on Zp c as open circles with arrows towards higher values. Iron abundances inferred by both studies are mostly consistent, 
and follow the 1:1 trend, albeit with significant scatter. Ath07 report systematically lower luminosities than HB06 for the lowest luminosity galaxies. Our 
conclusions are, however, robust to these uncertainties. 
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